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A simple method is presented to evaluate the effects of short-range correlations on the momentum 
distribution of nucleons in nuclear matter within the framework of the Green's function approach. 
The method provides a very efficient representation of the single-particle Green's function for a 
correlated system. The reliability of this method is established by comparing its results to those 
obtained in more elaborate calculations. The sensitivity of the momentum distribution on the 
nucleon-nucleon interaction and the nuclear density is studied. The momentum distributions of 
nucleons in finite nuclei are derived from those in nuclear matter using a local-density approximation. 
These results are compared to those obtained directly for light nuclei like le O. 
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I. INTRODUCTION 

Realistic nucleon-nucleon (NN) interactions like the phenomenological Reid soft-core potential jl| or One-Boson- 
Exchange (OBE) potentials ||, which are adjusted to fit the NN scattering data, typically contain rather strong 
short-range components. These short-range parts as well as a non-negligible tensor component are responsible for 
the fact that simple mean-field or Hartree-Fock (HF) calculations of nuclear systems yield very unsatisfactory results. 
It turns out that HF calculations using such realistic NN forces may not even lead to bound nuclei ||. Therefore, 
based on these theoretical considerations it seems obvious that nuclear wave- functions must contain correlations, 
which are induced by these short-range and tensor components and cannot be accounted for in the mean-field or HF 
approximation to the solution of the many-body problem. 

The question is, whether there exist experimental observables which reflect these correlations in an unambiguous 
way. In particular it would be nice if one could explore these correlations in terms of single-nuclcon observables 
since they are easier to measure as well as to calculate. This leads to the question how correlations affect the single- 
particle density p(r,r') in the nuclear many-body system. Rather than discussing this non-local representation of the 
density matrix, we may consider as well its Wigner transform [Q f{R, k). Integrating this Wigner distribution over all 
momenta k yields the local density p(r = r'). This local density distribution, or to be more precise the corresponding 
charge distribution has been investigated with high precision in elastic electron scattering experiments f| . Also the 
matter distribution can be analyzed in a rather model-independent way by means of elastic a scattering and other 
probes It seems, however, that these "experimental" charge and matter distributions can very well be reproduced 
within a mean-field approximation for the nuclear wave-function. 

Integrating the Wigner transform of the one-nucleon density matrix over all spatial coordinates one obtains the 
momentum distribution n(k). For an infinite system, invariant under local transformations, the mean-field or Hartree- 
Fock prediction for this momentum distribution is identical to the momentum distribution of a free Fermi gas. This 
means that all states with momenta less than the Fermi momentum kp are occupied with a probability n = 1, while 
all states with momenta k above kp are completely unoccupied [n = 0). Correlations beyond the HF approach modify 
this momentum distribution in the sense that states with momenta below kp are partly depleted, whereas states with 
high momenta are partly occupied. 

From these considerations for the infinite system of nuclear matter one may expect that correlations beyond HF 
will enhance the momentum distribution at high momenta k also for finite nuclei. This is one reason why modern 
electron accelerators have been used to explore the momentum distribution of nucleons in nuclei by means of nucleon 
knock-out, (e, e')p, experiments J?]-[| . 

Microscopic nuclear structure calculations which account for the effects of short-range and tensor correlations of 
realistic NN interactions are mainly performed for very light nuclei JlJ3|— p~3|] or infinite nuclear matter [^4|-pl|]. From 
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these results for nuclear matter one then tries to extract the effects of NN correlations in order to estimate their 
influence on the momentum distribution of real nuclei using a local density approximation (LDA) |^2|-^4|| . 

Recently, there have also been attempts to determine the momentum distribution in a microscopic calculation 
considering directly finite nuclei as 16 O |2^-|2^]. It turns out that the Green's function approach is particularly useful 
for these investigations. This method not only provides the total momentum distribution but also yields detailed 
information on the spectral function. This spectral function contains the information at which excitation of the residual 
nucleus or, using the nomenclature of the knock-out experiment, at which missing energy the various components 
of the momentum distribution should be observed. These studies predict that high-momentum components in the 
momentum distribution due to short-range correlations should show up preferentially at large missing energies. 

The comparison of the momentum distributions obtained within the Green's function approach for 16 O with cor- 
responding ones for nuclear matter exhibited remarkable discrepancies p7[ |. Therefore the question arises if these 
discrepancies demonstrate the limitation of the LDA in predicting the momentum distribution. As it has been argued 
already in p7[ | , such a conclusion would be premature since that comparison was plagued with various inconsistencies 
like e.g. 

• The momentum distribution calculated for nuclear matter using the Green's function approach has been available 
only for one specific density, the saturation density of nuclear matter. This density may be too large to be typical 
for the situation of nuclei as small as le O. 

• The result for nuclear matter has been derived from a self-consistent calculation of the single-particle Green's 
function Jl8| whereas the calculation for the finite system has been performed considering contributions to the 
self-energy of the nucleons up to second order in a nuclear matter G-matrix |2^ ] 

• The calculation of the self-energy for the nucleons in the finite system has been made employing a single-particle 
spectrum with a substantial gap at the Fermi surface, whereas a continuous prescription has been used for nuclear 
matter. 

• The calculation for the finite system was limited to partial waves with angular momenta I < 3. This limitation 
may be too severe for the momentum distribution at high momenta. 

• While the Reid soft-core potential (IJ has been used for the study of nuclear matter in []l8| the OBE potential 
B of 0] has been employed for the calculation of le O. 

It is one aim of the present investigation to remove some of the differences between these calculations of finite nuclei 
and infinite nuclear matter. Furthermore we want to study the sensitivity of the calculated momentum distributions 
in the nuclear systems on the nuclear density, the NN interaction considered and various other ingredients of the 
many-body calculation. For that purpose we have developed a new very efficient method to calculate the momentum 
distribution in nuclear matter using an approximation to the Green's function method very similar to the one presented 
in |27j for finite nuclei. The comparison of results obtained with this approximation scheme with those resulting from 
the much more sophisticated scheme of |l8| demonstrates the reliability of the method developed here. The efficiency 
of the new scheme allows the detailed studies mentioned above. 

After this introduction section 2 of this paper describes the technique to be used for studies of nuclear matter. In 
particular we will also present an efficient representation of the single-particle Green's function, which allows a self- 
consistent treatment. In section 3 we briefly review the basic approximations used in the calculation of the momentum 
distribution for finite nuclei of |p7|| and we outline a method to determine this momentum distribution, in which the 
mean field part is calculated for the finite system but the effects of correlations are taken from nuclear matter at 
various densities, using a LDA. The results of the numerical calculations are presented in section 4. In this chapter we 
discuss the sensitivity of the results in nuclear matter on the various ingredients. Also we compare the predictions of 
the LDA with results obtained by the method of [^7j . For that purpose we extended the studies of |27]] by considering 
different interactions and allowing for higher partial waves. The main conclusions are summarized in section 5. 



II. MOMENTUM DISTRIBUTION IN NUCLEAR MATTER 



A. Self-Energy and Dyson Equation 

Our calculation of the single-particle Green's function for nucleons in nuclear matter is based on the definition of 
the self-energy of the nuclcon, which includes the terms of first and second order in an effective interaction V, which 
we will define below. The expression for the term of first order, displayed in Fig Ilk, corresponds to the Hartree-Fock 
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expression for the single-particle energy of a nucleon with momentum k in a system of nuclear matter with a Fermi 
momentum k F 

^ HF \k)^yj <kh\V\kh> 

f ri\k-k F \ 
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In the second part of this equation the matrix elements of V are given using the conventional partial wave representation 
with L, S, J and T denoting the orbital angular momentum for the relative motion, the spin, the total angular 
momentum and the isospin of the two interacting nucleons, respectively. The relative momentum q is diagonal and 
an average value K av has been used for the cm. momentum, which is given by [28[ 

2 = f k 2 +q 2 , ifO<q< %\k-k F \ 

av \\k 2 -qk + \k 2 Fl ii\\k-k F \<q<\{k + k F ) {> 

The term of second order in the effective interaction with intermediate 2-particle 1-hole (2plh) states, displayed in 
Fig.reffig:diagb, also depends on the energy oj of the nucleon under consideration and can be calculated according to 
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(3) 



The single-particle energies e g correspond to the Hartree-Fock approximation for the single-particle energy 

«.-£ + *">(.), (4) 

with m for the mass of the nucleon. In the second part of eq.([|) we have used the so-called angle-averaged approxi- 
mation for the Pauli operator, which is defined by 128] 



0, if q<y/k*- 

Q(K,q)=i i K2 +£- k * , if y/k 2 F -KZ<q<k F +K (5) 

1, q iiq>k F + K 

The other contribution to the second order self-energy with intermediate 2-hole 1-particle states (2hlp), displayed in 
Fig.|l|c, can be calculated in a way very similar to eq.([|) 

yf yf <fcp|V|/n,/i 2 > 2 
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The Pauli operator Q, which ensured in eq.(||) that the sum over intermediate 2-particle states is restricted to states 
above the Fermi surface, is replace by a corresponding operator P to ensure that the intermediate 2-hole state are 
below the Fermi level 
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P(K,q) = 



0, 

k 2 F -K 2 -q 2 



1, 



2Kq 



if q>Vk 2 F - 

if - K 2 >q>k F -K 

if <7 < k F - K 



(7) 



The definitions of the mean value for the center of mass momentum K av and the weight function W(q) in the integral 
of eq.(^) are a bit more involved than in the case of the 2plh term and are given in the appendix | p9| . 

After the definition of the self-energy we could now proceed and calculate the corresponding single-particle Green's 
function by solving a Dyson equation of the form [p0| 



g(k, u) = g (HF ) (*,«)+ g (HF ) (k, u) [E< 2 * lh > (k, w) + (k, u)\ g{k, u) , 

with the single-particle Green's function in the Hartree-Fock approximation 

(HF)n.,.\_ ®{kp-k) t Q(k-k F ) 
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The momentum distribution n(k) can then be calculated from the imaginary part of the single-particle Green's 
function by 
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(10) 



B. Numerical Approach 



Instead of proceeding along the lines indicated in eqs.(Q) - (|1C|) , we use the fact that in all numerical calculations 
the integrals of eqs.(0) and (ra) will be discretized. This means that eq.(0) takes the form 
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while eq.(0) can be rewritten as 
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This discretization implies in particular that we represent the singularities of the self-energy in terms of discrete 
poles slightly above (2hlp) and below the real axes. This analytic structure of the self-energy is identical to the one 
obtained for a finite system within a model-space defined in terms of discrete single-particle states. This means that 
we may use the same techniques to determine the features of the single-particle Green's function as employed e.g. in 
pl|p2| . Translating this technique into the present example, this means that the single-particle Green's function will 
be defined in the Lehmann representation by 



N + M+l 



ui a ±ir] 
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with the sign in front of the infinitesimal imaginary part r] being positive for poles Lu a above the Fermi energy Ep and 
negative else. For each k the positions of these poles, oj a , and the residua, X 2 , can be determined from the solution 
of the following eigenvalue problem 
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Note that the dimension of this matrix (N + M + 1) as well as the matrix elements Fi, E^ 2plh ^ and Gj, E^ hlp ^ refer 
to the nomenclature employed in eqs.(^Tj) and (|l2|), respectively. Using the representation of the Green's function in 
eq.(|l^) the occupation probabilities are calculated easily as 

n(k) = ^e(E F -u a )Xl, (15) 

(X 

which leads to the momentum distribution if we divide by the single-particle density 
with kp the Fermi momentum of the nuclear matter system. 



C. BAGEL Approximation 



Using the Lehmann representation of the single-particle Green's function of eq. ( |13| ) the continua of states of nuclear 
matter with one additional nucleon and one hole are rep resented in terms of some discrete energies u a . Depending 
on the accuracy of the discretization on eqs.(jll|) and ( |12[ ) the number of eigenvalues typically considered in numerical 
calculations ranges from a few hundred up to a few thousand. This may be compared to the Hartree-Fock approxi- 
mation of eq.(^), in which the Green's function for a nucleon with momentum k is represented by just one pole. This 
number of pole terms is not a problem as long as one is just interested in the evaluation of the Green's function or 
simple observables as the momentum distribution. The structure of the Green's function in eq. (|l3]) , however, may 
be too complicate to be used in the evaluation of quantities which are defined in terms of products of these Green's 
functions. Examples for such quantities are e.g. the various response functions of nuclear matter or a self-consistent 
evaluation of the self-energy, which implies that the self-energies to be used in the Dyson equation (||) are calculated 
in terms of the resulting Green's functions. For such calculations it may be preferable to "optimize" the number of 
pole terms in eq.([H|), which means: Try to find a minimum number of poles, which yields the same observables than 
the complete Green's function. 

In order to develop such an efficient representation of the Green's function we try to apply the so-called "BAsis 
GEnerated by Lanczos" (BAGEL) scheme, which has successfully been used for the description of finite nuclei in 
finite model spaces [[53 3lJ 3^] . For that purpose we consider the operator a which corresponds to a part of the matrix 
in eq. (0) 
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and apply this operator on the single-particle state \a >, which in terms of the matrix representation of eq. ( |l7| ) 
described by the column vector (1,0... 0) T 



is 



a\a >— ek\a > +6i ai > , 



(18) 



where \a% > is orthogonal to \a > and the coefficient 5i is chosen so that \a± > is normalized. Following the 
Lanczos algorithm j34j], one can subsequently construct additional states \ai >, which are all orthogonal to each 
other. Applying the Lanczos procedure n times one obtains n basis states of the 2plh configuration space. 

In a similar way we can furthermore construct m basis states of the 2hlp configuration space by considering the 
corresponding sub matrix of eq.(14) 
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and reduce the eigenvalue problem of eq.(|l4j) to the corresponding one in the subspace defined by the basis of the 
single-particle state plus the (n + m) basis states generated by the Lanczos scheme just outlined. The Green's function 
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of this BAGEL (n,m) approximation is then defined according to eq.([13|) using the (n + m+1) eigenvalues and vectors 
obtained from the diagonalization of the matrix truncated to the subspace. It is obvious that the BAGEL(0,0) 
corresponds to the HF approximation, while for n approaching N and m close to M the BAGEL(n,m) approximation 
for the Green's function becomes identical to the exact solution of eqs. 



14) and @. 



D. Effective Interaction 



At the end of this section we want to define the effective NN interaction V, used in the definition of the self-energy 
above. One possible choice would be of course to replace V by the bare NN interaction. As it has been discussed 
already in the introduction, the HF approximation in terms of a realistic NN interaction is not a very useful approach 
and it is not clear whether a perturbation expansion in terms of the bare NN interaction up to second order, as 
just outlined, will be sufficient. Therefore we employ the G matrix, an appropriate solution of the Bethe-Goldstonc 
equation for V. The starting energy Z in the Bethe-Goldstone equation is chosen according to the Brueckncr-Hartree- 
Fock (BHF) choice for the self-energy of a nucleon with momentum k below the Fermi momentum and put to be the 
average of two single-particle states below the Fermi energy if k is above the Fermi momentum. This ensures that G 
remains real. With this choice we employ an approach which is very similar to the one used for finite nuclei in [g7| , 
where the self-energy is also calculated including terms up to second order in a nuclear matter G-matrix. 

Using the G matrix for the effective interaction also implies, however, that we have to face a double-counting 
problem. The diagram of second order, displayed in fig. [lp is to some extent already taken into account in the 
Brueckner-Hartree-Fock (BHF) approach for the self-energy displayed in fig.|l|a. This double-counting does not directly 
effect the calculation of the momentum distribution. The choice for the starting energy just presented leads to a real 
self-energy contribution of eq.([j]) without any poles and therefore the momentum distribution calculated for a self- 
energy, which only accounts for this term, remains identical to the HF one. There is a self-consistency problem, 
however, with respect to the energy spectrum reflected in the poles of the Green's function (see eq.(|l^)). If for the 
moment we ignore the 2hlp contribution to the self-energy and evaluate the Green's function according to scheme 
outlined in eqs.(|l3|) - ([IJ) for a nucleon with momentum below hp, we will find one eigenstate of eq.(|l4|) with 
negative energy and a large coefficient X a , the quasihole state, and N eigenvalues at positive energies. Due to the 
diagonalization, however, the energy of the quasihole state, e ^ h ' 2plh ^ w {\\ b e substantially below the corresponding 
HF energy e k . Therefore we replace the first element of the matrix in eq. (|l4|) by 

e fe => e k = e k - { t { ^ 2plh) - e k ) . (20) 

This shift in energy ensures that the quasihole state for a self-energy with inclusion of only the 2plh term will 
essentially be identical to the BHF energy. Therefore the double-counting is removed. Note again that this double- 
counting problem does not affect the calculation of the momentum distribution. The energy-shift is useful, however, 
to obtain a realistic energy spectrum for the poles of the Green's functions. 



III. MOMENTUM DISTRIBUTION IN FINITE NUCLEI 



A. Direct Approach 



directly for the nucleus 16 O is based on a 
self-energy of the nucleon calculated up to second order in a nuclear matter G-matrix as described by the diagrams 
of fig.|l]. As a first step one considers the HF contribution to the self-energy 

£^(fci,fci) = 2(2 1 ^ E (2J + 1)(2T+ 1) (hhjm^JTlGlkihhn^JT) . (21) 

^ Jl+ > n 2 l 2 j 2 JT 

The matrix elements of G used in this expression are antisymmetrized NN matrix elements calculated in the laboratory 
system The quantum numbers U and ji refer to the orbital and total angular momentum of the single nucleons in 
this frame and J and T denote the angular momentum and isospin of the 2-particle states. The matrix elements 
are calculated in a mixed representation with Uj referring to the radial quantum numbers of oscillator bound (hole) 
states, whereas the fcj denote the absolute value of the momentum for a free particle state. The summation over 
the oscillator quantum numbers is restricted to the states occupied in the independent particle model of 16 O. This 
Hartree-Fock part of the self-energy is real and does not depend on the energy. The HF single-particle wave functions 
can be obtained by expanding them 



Also the calculation of the momentum distribution presented in |E6 27 



G 



\a HF ljm >= \Kiljm >< K i \a ilF > tj (22) 

i 

in a complete and orthonormal set of regular basis functions within a spherical box of radius i?box which is large as 
compared to the radius of the nucleus 

$iOm(r) = (r\Kiljm) = N a ji{K t r)yi jm {$Lp) (23) 

In this equation 3^ym represent the spherical harmonics including the spin degrees of freedom and ji denote the 
spherical Bessel functions for the discrete momenta Ki which fulfill 



Using the normalization constants 



ji{KiR hox ) =0. (24) 



-, for/>0 



Na= v^-^*)' (25) 

^ -^L, for I = 0, 



the basis functions defined in Eq. (|23j ) are orthogonal and normalized within the box. The expansion coefficients of 
eq.(|2^) are obtained by diagonalizing the HF Hamiltonian 

£ m + E " F ( K ^ HF ) l3 = 4S (Ki\a H % • (26) 

n=l 

Here and in the following the set of basis states in the box has been truncated by assuming an appropriate -/V max . 
From the HF wave functions and energies one can construct the HF approximation to the single-particle Green's 
function in the box (compare eq.(@)), which for the finite nucleus has the form 

As an example for the contributions to the self-energy of second order in G we recall the calculation of the 2plh 
term. In the approach of p7j ] one first calculates the imaginary part of this self-energy contribution, depending on 
the energy to 

wffifaX;") = imrni En 2i2J2 J2 hUj3H £, T I Hdh J kjdk^j + i)( 2 t + 1) 

x (kihjiri2hj2JT\ G \k 3 l 3 j 3 k 4 l 4 j 4 JT) 
x {kzhhkdi34,J T \ GWxhh^khJT) 

x„6(u, + e^ a - (28) 

where the "experimental" single-particle energies e n2 i 2 j2 are used for the hole states (-47 MeV, -21.8 MeV, -15.7 MeV 
for S1/2, P3/2 arL d P1/2 states, respectively), while the energies of the particle states are given in terms of the kinetic 
energy only. The expression in Eq. (|2^) still ignores the requirement that the intermediate particle states must be 
orthogonal to the hole states, which are occupied for the nucleus under consideration. The techniques to incorporate 
the orthogonalization of the intermediate plane wave states to the occupied hole states as discussed in detail by 
Borromeo et al. |35|| have also been used here. The 2hlp contribution to the imaginary part Wf^ p {k\, k'^uj) can be 
calculated in a similar way (see also |35|). 

The choice to assume pure kinetic energies for the particle states in calculating the imaginary parts of W 2plh (Eq. 
(p8|)) and W 2hlp may not be very realistic for the excitation modes at low energy. Indeed a seizable imaginary part 
in W 2hlp is obtained only for energies uj below -40 MeV. As we are mainly interested, however, in the effects of 
short-range correlations, which lead to excitations of particle states with high momentum, the choice seems to be 
appropriate. A different approach would be required to treat the coupling to the very low-lying two-particle-one-hole 
and two-hole-one-particle states in an adequate way. Attempts at such a treatment can be found in Refs. |3(] |38[ ]. 

The real parts of the 2plh and 2hlp terms in the self-energy can be calculated from the corresponding imaginary 
parts by using dispersion relations p"4[ ]. As an example we present the dispersion relation for the 2plh part, which is 
given by 
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V h p h (k 1 ,k 1 ;w) = - 3 — — dw, (29) 



7T 



where P means a principal value integral. Putting the various contributions together the correction to the HF 
self-energy due to the second order terms can be written 

ASy (fa , fc 2 ; u) = (V 2plh -V c + V 2hlp ) + i (W 2plh + W 2hlp ) (30) 

where V c denotes a correction term to account for double-counting between the V 2plh and ladder contributions already 
contained in the HF part of the self-energy |27| (see also discussion at the end of section 2). With this correction to 
the self-energy one can solve a Dyson equation for the complete Green's function (see also eq.(^)), which corresponds 
to an integral equation for finite systems 



gi J (k 1 ,k 2 ;uj) = glf F \k 1 ,k 2 ;uj) + J dk 3 J dk 4 glf F \k 1 ,k 3 ;uj)Ai: lj (k 3 ,k 4 ;uj)gi :j (k 4 ,k2;uj). (31) 
From the imaginary part of this Green's function one can finally evaluate the momentum distribution according to 

n(fc)=^2(2j + l) f * dLu-Imag[ gij (k,k;Lo)] . (32) 



B. Local Density Approximation 

Instead of evaluating the momentum distribution directly for the finite nucleus one can try to deduce the effects of 
correlations on the momentum distribution from the investigation of nuclear matter. As a first step towards such an 
approach we consider the local density p HF (r) and the momentum distribution n HF (k) derived from the solution of 



the Hartree-Fock eq.(26). Using this density distribution we can define a local Fermi momentum 

37r 2 p HF (r) ll/3 



and evaluate an average occupation number for the states occupied in the mean field approach by 



(33) 



iV avcr = Jr 2 drp HF (r)n(k l ^ al (r)) , (34) 

with A the number of nuclcons (A=16 in our example of 16 0) and n the occupation number (see eq.(]l5|)) calculated 
for nuclear matter with the local Fermi momentum and averaged over all momenta below this Fermi momentum. 
With this average occupation number one can account for the depletion of the occupation of states occupied in HF 
approximation. The high momentum components originating from the partial occupation of states above the Fermi 
momentum are then evaluated as 

An(Jfe) = 4tt/c 2 J 47rr 2 drB(k~ k^ ca \r)) n (k; k l £ ca \r)) (35) 

where nik^kp] is the momentum distribution of nuclear matter according to eq.(|l6|) calculated for the local Fermi 
momentum. The total momentum distribution is then given as 

n LDA (k) = 7V avcr n HF (fc) + An(fc) . (36) 

Within this local density ap pro ximation we also would like to estimate the spectral strength which is missing in the 
calculation according to eq. (^|) due to a restriction in the sum of that equation to partial waves with orbital angular 
momentum up to I — L max . For our LDA approximation that restriction would mean to consider contributions to 
An(k) in eq.(B3) with 



\l\ = rk sin tp < TiyJ L max (L max + 1) , (37) 

ip denoting the angle between r and the momentum k. This means that the integrand in eq.(^|) should be reduced 
by a factor 



2 . Lmax\L"max ~t~ 1) /on\ 

— arcsm — . (38) 

7r rk 
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IV. RESULTS AND DISCUSSION 



A. Nuclear Matter 



In order to evaluate the momentum distribution for nuclear matter at a given density p, which may as well be 
characterized by the corresponding Fermi momentum kp, using the method outlined in the previous section, we 
have to determine as a first step the spectrum of HF single-particle energies e&. In our approach, defining the 
HF contribution to the self-energy (see [l]b) in terms of the nuclear matter G-matrix, these single-particle energies 
correspond to the single-particle energies obtained in the BHF approximation. For the OBE potential B of the 
BHF single-particle energies have been parameterized [p9[ in terms of an effective mass m* and a constant shift C by 



\Jk 2 + m* 2 -m* + m + C . (39) 



The parameters m* and C as a function of the Fermi momentum kp are listed in table 2 of [ p9[ . It should be mentioned 
that in our study we have used the so-called non-relativistic parameterization since in our present study we ignore 
all effects of the Dirac BHF approach due to a change of the Dirac spinors of the nucleons in the nuclear medium. 
With this definition of the single-particle spectrum one can evaluate the matrix in eq.(|l4|) with the renormalization 
of eq.(pOf), solve the eigenvalue problem of eq.(|l4|) and determine the momentum distribution with eq. (|l5|) . 

Results for this momentum distribution in nuclear matter at the empirical saturation density (/cp=1.36 fm _1 ) are 
displayed in Fig.|2|. The momentum distribution derived from the OBE potential (dashed-dotted line) is compared to 
the one obtained using the Reid soft-core potential employing the same technique (dashed line). The Reid soft-core 
potential predicts stronger effects of correlations in this momentum distribution. This is characterized by a stronger 
depletion of the states with momenta k below the Fermi momentum (the Reid soft-core yields an average occupation 
of these states of 0.83 while the OBE predicts 0.86 ) as well as larger probability at higher momenta (for k w 4 fm _1 
the density which is obtained for the Reid potential is by more than a factor 2 larger as the one deduced from the 
OBE). This is in agreement with the observation that the modern OBE potentials are "softer" and contain a weaker 
tensor force as it is also reflected in the D-state probability calculated for the deuteron 

Figure ^ also shows the prediction for the momentum distribution of nuclear matter obtained for the Reid soft-core 
potential using the much more sophisticated techniques of The good agreement of our approach, in which the 
self-energy is calculated in a perturbative scheme, including terms up to second order in G, with the one of |]l8| where 
the particle-particle hole-hole ladders are taken into account to all order using a self-consistent single-particle Greens 
function, gives us some confidence that the present approach provides reliable information for systematic studies in 
nuclear matter as well as finite nuclei. 

The sensitivity of the calculated momentum distribution on the nuclear density is demonstrated in fig||. In order 
to allow a direct comparison this figure does not show the momentum distribution but the occupation of the single- 
particle states as a function of the momentum k in units of the Fermi momentum (k/kp). The results for the 
occupation of states below the Fermi momentum (k/kp < 1) are shown with respect to the linear scale on the left 
axis, while the occupation of states with momenta larger than kp are shown with respect to the logarithmic scale on 
the axis at the right hand side. One observes that the calculated occupations are rather insensitive to density of the 
nuclear system. Only for very small densities (fci?=0.8 fm _1 , which corresponds to roughly 20 percent of the empirical 
saturation density) one finds occupation probabilities which are considerably smaller. This might be an indication of 
the instability of homogeneous nuclear matter at such small densities p3| , p4[ . 

In order to explore the density dependence a bit more in detail, we have separated the contributions to the single- 
particle density into a quasihole contribution and a continuum contribution. For momenta k below the Fermi momen- 
tum a large contribution to the momentum distribution of eq.(|l5|) originates from one eigenstate a with a maximal 
coefficient X a and an eigenvalue uj a which we identify as the quasihole energy E q h(k). In particular at small momenta 
(k sa 0.2kp) one also finds that a few states around the quasihole energy exhibit large coefficients X a and therefore 
contribute significantly to the sum in eq. ( |T5| ) . We define the quasihole strength to be the contribution of all terms in 
eq.(|l^) which originate from an eigenstate of eq.(14) with an eigenvalue within an interval of length 3 MeV around 



the quasihole energy E q h 

N qh (k) =J2 e ( E F~ ^a) &(E qh (k) + 1 .5 - w a ) e(u a - E qh (k) + 1.5) X 2 a . (40) 

a 

For momenta k larger than the Fermi momentum the eigenvalue of the state with maximal expansion coefficient 
X a occurs at energies above the Fermi energy Ep and therefore we don not obtain any quasihole strength for those 
momenta. The quasihole strength N q h is shown in fig.^ as a function of the ratio k/kp for two densities (lines labeled 







with triangles). The remaining contributions to the sum in eq.(|15J) will be called the continuum contribution to the 
occupation probability or momentum distribution. 

From the inspection of the results displayed in figj^ one finds that the quasihole contribution to the occupation 
probability increases drastically with the momentum while the continuum contribution decreases in a corresponding 
way. A typical ratio of quasihole versus continuum contribution is 0.6 for small momenta but as large as 10 for 
momenta close to the Fermi momentum. This means that the energy distribution of the single-particle strength 
is highly localized at the quasihole energy for states with momenta close to kp whereas one observes a brought 
distribution and fragmentation of the strength for small momenta, i.e. deeply bound hole states. 

The continuum contribution to the occupation probability decreases monotonically with increasing momentum and 
is a rather smooth function even at the Fermi momentum. This implies that the gap in the momentum distribution 
exhibited e.g. in fig|2] at k = kp originates simply from the fact that the quasihole contribution vanishes since the 
energy of the corresponding state gets larger than the Fermi energy. It is worth noting that about 70 percent of the 
single-particle strength is located in the quasihole contribution. Two thirds of the remaining continuum contribution 
occurs at momenta below the Fermi momentum and only one third of the continuum contribution, which means 
slightly more than 10 percent of the total strength occurs at momenta above kp. 

From the single-particle Green's function we also determine the mean value of the energy for the spectral distribution 
at a given momentum 

— ( tF duj wlmag qik.ui) 

= Mk) ' (41) 

or translated into the tools we are using in our numerical treatment 

Ct 

with the occupation number n calculated following eq(|l5|). The summation in this equation can be truncated as 
discussed above to determine the mean value for the energy resulting from the continuum part of the momentum 
distribution. 

Such mean values are presented in the right part of fig.^[ One finds that these mean values for the continuum are 
more negative than the corresponding HF single-particle energies or the energies of the quasihole states, which are 
shown in the left part of fig.^j (note the different scales on the axes). Particularly at large momenta, above kp where 
the continuum part represents the total momentum distribution, these mean values oj(k) are very attractive. This 
implies that these high momentum components of the momentum distribution occur predominantly at large excitation 
energies of the residual nuclear system, which corresponds to large missing energies in knock out experiments. This 
result for nuclear matter confirms the observations made for finite systems in |p6| , p7j| . 

In the left part of fig|| the BHF single-particle energies (lines labeled with triangles) are compared to the energies 
E q h of the quasihole energies. One can see that the inclusion of the 2hlp terms of fig.|l|c yields a repulsive contribution 
to the quasihole energy. This is specially true for states with momenta well below the Fermi momentum. This means 
that the removal energy for nucleon knock out experiments exciting states with large spectroscopic factor should be 
much smaller for these deeply bound states than predicted in BHF calculations. 

As a last point in this subsection we would like to explore the sensitivity of the calculated momentum distribution 
on the HF single-particle spectrum. For that purpose we have modified the effective mass parameter m* in the 
parameterization of eq.(||) from the BHF value of m*=623 MeV [|9) at k F = 1.35 fin" 1 by ± 100 MeV. From the 
upper part of fig^ one can see that a reduction of the effective mass, i.e. the single-particle energy shows a stronger 
momentum dependence, yields a reduction of the correlation effect. This reduction of the correlations is indicated 
by an enhancement of the occupation of states with k < kp and a reduction of the high momentum components. 
Another modification of the HF single-particle spectrum can be obtained by introducing a gap between the energies of 
particle and hole states. As on can observe from the lower part of figH such a gap does not affect the high momentum 
components very much but only the occupation probability around kp. 



B. BAGEL Approximation 



In order to investigate the efficiency of the BAGEL approximation introduced in the previous section for the 
representation of the single-particle Green's function in terms of a few poles we have evaluated the momentum 
distribution of nuclear matter considering various combinations BAGEL (n,m) for the number of basis states n in the 
2plh and m in the 2hlp part of the eigenvalue problem eq.(p^). Results for a few examples are displayed in fig|7]. 
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One finds that a very good approximation for the high-momentum components is obtained already with a very 
small number of basis states. The occupation probabilities for k > kp are reproduced in quite a satisfactory way in 
the simplest approximation BAGEL (1,1) and the results become indistinguishable from the exact results if we use 
any approximation with n, m larger than 1. The convergence of the BAGEL approximation towards the exact result 
with increasing n, m is not as good for the occupation of states below kp. A larger number of pole terms is required 
in particular to reproduce the decrease of the occupation number with k getting close to kp. 

Of course it is still a very efficient approximation to reduce the number of poles in the Lehmann representation 
of the single-particle Green's function of eq. ( |l3| ) from a few hundred obtained by an optimized discretization of the 
integrals in eqs.(|^) and (||) to n + m + 1 = 21 in the BAGEL(15,5) approximation, but a closer inspection may help 
us to reduce the number of terms even more. 

Analyzing the basis states which are generated by the BAGEL approach in the 2plh part of the Hilbert space 
by applying a of eq.(|l7]) according to the Lanczos scheme, one observes that states are generated with very large 
eigenvalue uj a but negligible amplitude X a . It requires some iteration steps to generate a few states with lower energy 
and non- negligible coefficient. This generation of basis states with extreme energies is of course a feature of the 
Lanczos approach, which is not optimal for our present purpose. The situation is better for the generation of basis 
states in the 2hlp sector, as the eigenvalues of A (see eq. 17|)) are more limited. 

From this discussion we see that the BAGEL approximation for the Green's function can be made more efficient 
by either ignoring the contributions of those poles in the Lehmann representation of eq.([l3]), which show very small 
coefficients X a or to replace the Lanczos algorithm of generating the basis states by one, which preferably generates 
eigenstates close to the Fermi energy. 



C. Finite Nuclei 



As a first example for the momentum distributi on cal culated directly for the finite nucleus 16 O using the method 
described in p7f and briefly reviewed in section III A , we present in fig.|| as a typical example the momentum 
distribution obtained for Pin partial wave. This refers to the corresponding contribution to the sum in cq.(|32]) 
without the factor 2(2j + 1) for degeneracy of these states. As in J27J we split the momentum distribution into a 
quasihole contribution, which should be observed if the {A — 1) nucleus remains in its groundstate, and a continuum 
part reflecting the momentum distribution observed at larger missing energies. 

The results obtained for the Reid soft-core potential pj (dashed lines) are rather similar to those evaluated for the 
OBE potential B of Q (solid lines). However, there are some characteristic differences which can also be observed in 
the other partial waves: The quasihole contribution evaluated for the Reid potential exhibits a maximum at smaller 
momenta and drops faster with increasing momentum. This reflects the fact that nuclear structure calculations like 
BHF yield less binding energy and a larger radius using the Reid soft-core potential as compared to the OBE model 
for the NN interaction. The continuum part, on the other side, exhibits larger contributions at high momenta using 
the Reid potential. Following the arguments present in the introduction of this paper, this would be an indication 
that the Reid potential predicts "stronger" correlations. 

Another indicator for the importance of NN correlations are the occupation probabilities for the various partial 
waves as they are listed in table |l]. These occupation probabilities are obtained by a momentum integration of the 
various partial wave contributions in eq.(|3^). Again we distinguish between quasihole and continuum contribution 
for the partial waves with I < 1 and compare the results for the two models of the NN interaction. The results 
presented here for the OBE potential deviate slightly from those presented in |27|] as we have increased the interval 
for the energy integration in eq.(|3^) in our present study. It is interesting to note that the occupation probabilities 
are quite similar for both interactions for the orbits with angular momentum I < 2 with slightly larger values for 
the OBE potential. For the partial waves with larger I, however, the Reid potential predicts occupation probabilities 
which are significantly larger than those for the OBE potential. This difference seems to be due to the stronger tensor 
component contained in the Reid potential. 

Multiplying the occupation probabilities of table | with the degeneracy factors 2(2j + 1) one finds that 2.05 (2.09, 
using Reid) "nucleons" out of the 16 for 16 are represented by the continuum part of the momentum distribution. The 
total nucleon numbers, including the quasihole part are 16.07 and 15.96 for the OBE and Reid potential, respectively. 
This means that the particle-number violating features of our present approach lacking a self-consistent treatment of 
the single-particle Green's function are not very strong [|o) . 

The continuum part of the total momentum distribution including partial waves with / < 4 in the summation of 
eq.(|3^) is displayed in the left part of fig. ||. Again we can observe the characteristic differences obtained for the two 
interactions: While the OBE potential yields a momentum distribution which is slightly larger at small momenta, the 
Reid potential predicts contributions which are larger by almost a factor of 2 at large momenta. This is in complete 
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agreement with the observation made above for nuclear matter (see fig]|). 

Figure ^ also exhibits the spectral distribution of the continuum part. This spectral distribution is obtained by 
replacing the energy integration in eq.(32) by a momentum integration. One observes that the spectral distribution 
derived for the OBE potential is slightly larger at energies with small absolute value, which corresponds to small missing 
energies, while the Reid potential yields a larger result at large missing energies. Concluding we may characterize 
the differences between the two interactions by the statement: The Reid potential predicts a larger component of the 
single-particle density at large momenta and large missing energies as compared to OBEP B. 

Finally, we would like to discuss the validity of the local density approximation (LDA) introduced in section [II B for 
the description of the momentum distribution at large momenta. Since this high-momentum part of the momentum 
distribution is dominated by the continuum contribution An(fc) in eq.(|36|) we will restrict our discussion to this part 
and its comparison to the continuum part of the momentum distribution evaluated directly for finite nuclei. 

The result for An(k) calculated according to eq.(|35|) using the OBE interaction is represented by the solid line in 
the left part of figjlO[ In the discussion of the LDA above we already introduced a scheme to simulate the effects of 
a cutoff in the partial wave summation for the evaluation of the momentum distribution in finite nuclei. Comparing 
the prediction for the total An(fc) with those employing a restriction to a maximal orbital angular momentum, one 
observes that such a restriction leads to remarkable differences in particular at high momenta. Using L max — 3 in 
eq.(|3^) one obtains a prediction for the momentum distribution, which is only one half of the result including all 
partial waves at k « 3.5 fin . 

The right part of fig.|l^ shows the comparison of the continuum part of the momentum distribution evaluated 
directly for 16 O with the LDA prediction using L max = 4, the maximal orbital angular momentum, which we have 
taken into account in our direct evaluation for the finite system. The agreement of the LDA with the direct evaluation 
is very good in particular at high momenta. For small momenta the direct calculation for the finite system yields 
a larger value than the LDA. For these small momenta, however, one must keep in mind that the continuum part 
calculated directly for the finite system contains components, which have a momentum distribution similar to the 
quasihole part 27|, whereas the An(fc) evaluated in LDA only accounts for momenta above the local Fermi momenta 

MO- 

Concluding we may say that the discrepancy observed in ref. pTfl between the predictions of high momentum 
components in the single-particle density of nuclear matter and those for finite systems has been resolved. The LDA 
seems to produce very similar results if the same NN interaction is used, an appropriate average over nuclear matter 
with various densities is considered and the effects of truncating the partial wave expansion in finite systems are taken 
into account. 



V. CONCLUSIONS 



The single-particle momentum distribution has been investigated for nuclear matter and the finite nucleus employing 
two different realistic models for the NN interaction. The investigations are based on the Green's function approach 
approximating the self-energy of the nucleon including all contributions up to second order in the G-matrix. The 
main results can be summarized by the following conclusions. 

• The present approach yields results for nuclear matter which are in very good agreement with the more elaborate 
calculations of [ fl8|| . This gives us some confidence that the same approximation should also produce reliable 
results if applied for finite systems. 

• The momentum distribution observed in knock-out experiments with small missing energies should mainly 
observe the quasihole distribution with small components at high momenta. Larger contributions to the single- 
particle density at high momenta should be obtained at large excitation energies of the residual nucleus. 

• The prediction for the momentum distribution depend weakly but in a rather characteristic way on the inter- 
action used. The stronger tensor and short-range components of the Reid soft-core potential yields a larger 
single-particle strength at high momenta and large missing energies as compared to the OBE potential B of 

. For momenta around 3.5 fm the momentum distribution derived from the Reid interaction is larger by a 
factor 2 in nuclear matter as well as finite nuclei. 

• The momentum distribution calculated for finite nuclei is rather sensitive to a truncation in a partial wave 
expansion. Orbits with angular momenta I larger than 4 should be taken into account to obtain stable results 
at momenta k ps 4 fm . 
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• A local density approximation, in which the high momentum components in the single-particle density are 
derived from the study of correlations in nuclear matter, yields a very good agreement with corresponding 
studies for finite nuclei, if the effects due to truncations in the partial wave expansion are considered. 

• The numerical scheme developed for the solution of the Dyson equation in nuclear matter leads to a very efficient 
representation of the single-particle Green's function in terms of a few "characteristic" poles in the Lehmann 
representation. This BAGEL approximation could be very useful e.g. in studies of nuclear response functions 
beyond the HF and RPA approximations. 
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hospitality at the Facultat de Fisica, Universitat de Barcelona, and the support by the program for Visiting Professors 
of this university. 



VI. APPENDIX 



This appendix lists the expressions p9[ for the weighting function W(q) and the average cm. momentum K av 
which have been used in the calculation of the 2hlp contribution to the nucleon self-energy according to eq.(^). For 
that purpose we distinguish 4 different cases of the momentum k for which the self-energy shall be evaluated: 

(a) < k < \k F 



case (a) 



{b)\k F <k <k F 

(c) k F <k < 3k F 

(d) 3k F < k < oo . 



' !((fc + ?) 3 + (!(* 2 + *!)-9 2 ) § ) 

\{k 2 -k 2 F )+q{q + k) 

K av (q,k) = \ H(g + fc ) 3 -(g- fc ) 3 ) 

2kq 

i((M 3 -(<z-£0 3 ) 



\{k 2 F -k 2 -q 2 ) + kq 



for Uk F -k) <q< \{k + k F ) 



for \{k + k F ) <q<(k F -k) 
for (k F - k) < q < (k + kp) 



(43) 



(44) 



case (b) 



Q ■ {\(k 2 - k%) + q(q + k)} fori(fc F - k) < q < \{k F + k) 
W(q) = \ q ■ {2kq} for±(/fc F + k) < q < (k F - k) 

q ■ {i(k 2 F -k 2 - q 2 ) + kq} for(fc F - k) < q < {k F + k) 



' 1 h k + q f + { l {k 2 + k 2 F) _ q 2 )i \ 

V " L for \{k F - k) < q < (k F - k) 



K av (k,q) = < 



\(k 2 -k 2 F ) + q(q + k) 

\(kl-{\{k 2 + k 2 F )-q 2 )^) 
\(k 2 p - k 2 ) + ^q 2 

l((fc F ) 3 -(g-fc) 3 ) 
[ \{k 2 F -k 2 -q 2 )+kq 



for (k F - k) < q < \{k + k F ) 
for Uk + kp) < q < (k + kp) 



(45) 



(46) 
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' q ■ {j(k 2 - k 2 F ) + q(q + k)} for±(fc F - k) < q < (k F - k) 
W(q) = I k ■ {\{k 2 F - k 2 ) + \q 2 } for(fc F - k) < q < \{k F + k) 

k k ■ {\{k 2 F -k 2 - q 2 ) + kq} for±(fc F + k) < q < (k F + k) 



case (c) 



' \(k%-{\{k 2 + k 2 F )- q 2 )i) 



K av {k,q) = < 



I(fc|-fc2 )+ l 9 2 

i((M 3 -k-fc| 3 ) 

I \{k 2 F -k 2 -q 2 ) + kq 



for yj±(k - k F )(k + k F ) <q< \{k + k F ) 
for \{k + k F ) <q<(k + k F ) 



W(q) = { 



' q ■ {j(k 2 F - k 2 ) + ±q 2 } fory^fc - k F )(k + k F ) < q < \{k F + k) 
k q ■ {\{k 2 F - k 2 - q 2 ) + kq} for±(/c F + k) < q < (k F + k) 



case (d) 



Kav(k,q) 



i((M 3 -|g-fc| 3 ) 

\{k 2 F -k 2 -q 2 ) + kq 



for (k~k F ) <q< (k + k F ) 



W(q) = { q ■ {\{k% - k 2 - q 2 ) + kq} for(fc - k F ) < q < (k F + k) 



(47) 



(48) 



(49) 
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TABLE I. Distribution of nucleons in 16 0. Listed are the occupation probabilities for various partial waves. These prob- 
abilities are obtained by integrating the partial wave contributions in eq.(32) ignoring the degeneracy factors 2(2j + 1). For 
states with / < 1 the contributions from the quasihole (n qh ) and the continuum part (n c ) of the spectral function are listed 
separately. Results have been obtained using the OBE potential B and the Reid soft-core potential. 



Ij 




OBEP B 


n c 




Reid 


n c 


Sl/2 


0.780 




0.157 


0.778 




0.117 


P3/2 


0.914 




0.042 


0.896 




0.040 


Pl/2 


0.898 




0.046 


0.896 




0.047 


d-5/2 






0.022 






0.018 


^3/2 






0.027 






0.024 


/7/2 






0.008 






0.014 


fa/2 






0.013 






0.019 


59/2 






0.002 






0.004 


57/2 






0.004 






0.007 
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FIG. 1. Graphical representation of the Hartree-Fock (a), the 2-particle 1-hole (2plh, b) and the 2- hole 1-particle contribution 
(2hlp, c) to the self-energy of the nucleon 

FIG. 2. Momentum distribution in nuclear matter at the empirical saturation density, &f =1.36 fm _1 . Results obtained 
with the approximations discussed in section 2 for the OBE potential B (dashed-dotted line) and the Reid soft-core potential 
(dashed line) are compared to the results using the Reid potential (solid line "exact") as derived from the more sophisticated 
calculations of ref.[18] 

FIG. 3. Occupation probabilities in nuclear matter as a function of the momentum in units of the Fermi momentum &f- 
Results obtained for the OBE potential B are shown for various densities (see description in the figure). The occupation 
probabilities for momenta below &f are displayed with reference to the axis on the left side of the figure, while those for k 
larger &f refer to the logarithmic scale on the axis at the right side. 

FIG. 4. The contribution of the quasihole state to the occupation probability (curves labeled with triangle, see eq.(40)) and 
the continuum contribution as a function of the momentum k/kF calculated for the OBE potential B. The dashed lines refer 
to a Fermi momentum of nuclear matter of 1.0 fm _1 , whereas the solid lines are obtained for &f = 1-35 fin -1 . 

FIG. 5. Energy spectra for nuclear matter with Fermi momentum fcf = 1 fm _1 (dashed curves) and &f = 1.35 fm _1 (solid 
lines) as a function of momentum. The left part of the figure exhibits the HF single-particle energies (curves labeled with 
triangles) and E q h the energies of the quasihole states. In the right part the mean values uj of the continuum part of the 
momentum distribution (see eq.42) are displayed. 

FIG. 6. Occupation probabilities in nuclear matter at /cf=1.35 fm _1 using various modifications of the BHF single-particle 
spectrum 

FIG. 7. Occupation probabilities in nuclear matter at /cf=1.35 fm _1 calculated in various BAGEL(n,m) approximations are 
compared to the result obtained with the complete Lehmann representation of the single-particle Green's function. 

FIG. 8. Momentum distribution for O in the P\/2 partial wave (see eq.(32)). The distributions are normalized such that 
/ dk n(k) = 1 if one orbit would be occupied. The three parts of the figure display the quasihole contribution, the continuum 
contribution and the sum of these two (total) as obtained for the OBE B (solid lines) and the Reid soft-core potential (dashed 
lines). 



FIG. 9. Continuum part of the single-particle density for m O as a function of the momentum (left part, see eq.(32)) and 
energy (right part). Results are presented as obtained for the OBE B (solid lines) and the Reid soft-core potential (dashed 
lines). 

FIG. 10. Continuum part of momentum distribution for 1(i O as obtained from the local density approximation (see eq.(35)) 
using the OBE potential B. In the left part of the figure the total contribution is compared to results obtained assuming 
various values L max for the truncation of the partial wave expansion. The right part of the figure compares the local density 
approximation with L max —4: with the result obtained in a direct calculation of 1(i O. 
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